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Abstract 

Recent variational calculations of the deuteron and the triton illustrate that simple 
wave function ansatze become more effective after evolving the nucleon-nucleon 
potential to lower momentum ("VjowA:")- However, wave function artifacts from the 
use of sharp cutoffs in relative momentum decrease effectiveness for small cutoffs 
(< 2fm~^) and slow down convergence in harmonic oscillator bases. These sharp 
cutoff artifacts are eliminated when Fiowfc is generated using a sufficiently smooth 
cutoff regulator. 



1 Introduction 



Variational calculations for nuclei are complicated by the strong short-range 
repulsion and tensor forces of conventional nucleon-nucleon potentials, which 
necessitate highly correlated trial wave functions. However, the nonpertur- 
bative nature of conventional inter-nucleon interactions is strongly scale or 
resolution dependent and can be radically modified by using the renormaliza- 
tion group (RG) to lower the momentum cutoff of the two-nucleon potential 
[1]. A particular consequence is that the short-range correlations in few- and 
many-body wave functions are significantly reduced at lower resolutions [1]. 
This has the practical implication that variational calculations should be more 
effective at lower cutoffs using simple wave function ansatze. A recent study us- 
ing low-momentum potentials (generically called "Viowfc") supports this claim 
[2]. This conclusion is also consistent with the results of Viviani et al. [3], 
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who showed that variational calculations of 3- and 4-body nuclei using hy- 
perspherical harmonics converge faster for potentials with greater nonlocality, 
which reduces short-range correlations. 

The optimistic conclusions regarding the use of VJowfe in variational calcula- 
tions are clouded by some problems associated with the sharp momentum 
cutoff, which are of concern for some practical applications. In particular, one 
would expect that for very- low-energy observablcs such as the deutcron or 
triton binding energies, there should be improved variational estimates down 
to very low cutoffs (e.g., comparable to the "binding momentum"), and rapid 
convergence with basis size. However, what was observed in Ref. [2] was im- 
provement down to moderate cutoffs of 2 fm^^ followed by a degradation of 
the variational estimates at smaller cutoffs. Further, for a wide range of cut- 
offs the convergence of the triton binding energy with the size of a harmonic 
oscillator basis was exceedingly slow once the energy prediction was at the 
100 keV level. Both of these problems were attributed to the use of a sharp 
cutoff on the relative momentum. In this letter, we revisit the variational cal- 
culations in Ref. [2] to demonstrate that these limitations can be eliminated 
by constructing Viowfc potentials using a smooth regulator. There are many 
new issues to consider with such a regulator, but we defer most of the further 
discussion to a more complete and wide-ranging investigation [4] . 



2 Mowfe with a Smooth Regulator 



The construction of VJowfc with a sharp cutoff is documented in Refs. [5,6], 
where it is shown that either RG equations or Lee-Suzuki transformations can 
be used. The latter approach relies heavily on the introduction of orthogonal 
low- and high-energy subspaccs with projection operators P and Q, such that 
P + Q = 1 and PQ = QP = 0. In momentum space for the two-nucleon 
system, the last condition mandates a sharp cutoff A in relative momentum, 
so that P-space integrals run from to A while Q-space integrals run from 
A to oo (or to a large bare cutoff). But while replacing a sharp cutoff with a 
regulator that smoothly cuts off the relative momentum seems incompatible 
with methods requiring PQ = 0, it is not a conceptual problem for the more 
general RG approach [5]. Details will be presented in Ref. [4]; here we simply 
summarize a three-step procedure applied in the present calculations. 

Smooth cutoff regulators will be applied in each partial wave as simple func- 
tions of the relative momentum. It is convenient and efficient for numerical 
calculations to define the smoothly regulated energy-independent potential 
Viowk s-nd the corresponding Tiowfc matrix in terms of a reduced potential v 
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and a reduced T matrix t as 

Vio^k{k',k) = f{k')v{k',k)f{k) , (1) 

and 

Tiowfe(A:', k; = f{k')t{k\ k; k')f{k) , (2) 
where f{k) is a smooth cutoff function. Here we adopt for /(/c) the exponential 
form used in chiral EFT potentials at N^LO order [7], 

f{k) = e-^'^'/^')" . (3) 

We use n — A throughout this work. The reduced t matrix obeys a Lipmann- 
Schwinger equation with loop integrals smoothly cut off by the internal factors 
of/(p) [5], 

t{k ,k;k ) = v{k ,k) + - p dp . 4) 

n Jo K'^ — p^ 

Note that the cutoff is on the loop momentum but not on external momenta. 
Principal value integrals are implicit throughout. 

In the energy-independent RG approach, an RG equation for the reduced 
interaction v{k', k) is derived by demanding that -^t{k\ k; k"^) = 0. ^ Using the 
large-cutoff initial condition v{k',k) — V^Nik'^k), the resulting set of coupled 
differential equations can be numerically integrated to evolve the interaction 
to lower cutoffs. The resulting Viowfe preserves the original on-shell T matrix 
up to factors of the smooth regulator, 

Txo^kiK k; k^) = f\k)T^N{K k; k^), (5) 

which implies that low-energy phase shifts are preserved for on-shell momenta 
away from the transition region near the cutoff where the regulator function 
rapidly decreases to zero. While the energy-independent RG approach provides 
a direct path to construct the smooth cutoff version of VJowfc, it is not the most 
convenient or numerically robust method [4]. 

In this work, we prefer to use a much simpler energy- dependent RG equation 
to evolve the bare potential Vnn to a lower cutoff. The energy-dependent RG 
equation is obtained by requiring invariance of the full off-shell T matrix, 
■^t{k', k; E) = 0, which can be formally integrated [4] to recover the Bloch- 
Horowitz equation with a smooth cutoff: 

, , 2 p 2 , (. f2( ^\ VNN{k',p)v{p,k;E) 
v{k , k] E) ^ Vf^^{k , k) + - pdp\l-f{p)\ . (6) 

TT Jo \ I tj — p 

^ This RG equation generates a non-hermitian interaction. Hermiticity can be re- 
stored by using a symmetrized form of the RG equation [4] . The resulting interaction 
preserves the on-shell T matrix, but no longer preserves the half-ofFshell T matrix. 
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The resulting integral equation for v{k', k; E) is much more efficient for nu- 
merical calculations than the original set of coupled differential equations. 

The second step is to trade the energy dependence for momentum dependence 
by defining an energy- independent (but non-hermitian) Viowfe(^', k) that gives 
the same half-on-shell T matrix and wave functions as the energy-dependent 
interaction V^sik', k; E) = f{k')v{k', k; E)f{k), 

(fc'|Kfr(/)|^p) = (A^-'|Mowfc|^p) . (7) 

The are the self-consistent wave functions of the energy-dependent low- 
momentum Hamilitonian Hcs{p^)- Using the completeness of these wave func- 
tions, one obtains a simple expression for the non-hermitian Viowfc, 

-\ I p'dpj^ k"'dk"V,s{k',k";p')%{k")%*{k) , (8) 

where ^* is the biorthogonal complement wave function. Note that the integral 
over the continuous scattering states will include a summation over discrete 
bound states, when present. 

The final step is to apply a Gram-Schmidt procedure to hcrmitizc the poten- 
tial, as prescribed in Ref. [8]. The end result is a hcrmitian, energy-independent 
Viowfe with a smooth cutoff regulator that preserves the low energy on-shell T 
matrix up to factors of the regulator as in Eq. (5). In addition, the complete 
set of wave functions obtained from diagonalizing Hio^ k can be used to con- 
sistently evolve general operators beyond the Hamiltonian with the smooth 
cutoff. Thus, the present approach may be viewed as a generalization of con- 
ventional effective interaction methods such as Lee-Suzuki transformations to 
smooth cutoffs. Further details concerning the present approach as well as the 
energy-independent RG method will be provided in Ref. [4]. 

As in Ref. [2], we will show results starting from the Argonnc ?;i8 potential [10] 
as Vnn, since it has been used in almost all modern calculations of light nu- 
clei. However, as with the sharp cutoff Viowfc calculations, the pattern of results 
for the full cutoff range shown here does not vary significantly with different 
initial potentials. When using the sharp cutoff Viowfc, two-body bound-state 
properties and phaseshifts arc preserved by construction for external relative 
momenta right up to the cutoff. (Note that three- and many-body observ- 
ables require the consistent addition of a three-body force to remove cutoff 
dependence [9].) The use of a smooth regulator, however, distorts the phase- 
shifts near the cutoff according to Eq. (5) to a degree that depends on the 
type of regulator function [4]. These distortions are not important for the low- 
energy observables discussed here, but will need further assessment for future 
applications. 
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Since the final Viow k is energy independent and hermitian, variational calcula- 
tions proceed as described in ordinary quantum mechanics texts (e.g., without 
special normalizations as needed for energy-dependent potentials). That is, 
given a trial wave function '0triab our variational estimate for the ground state 
energy at cutoff A is: 

p /AN (V^trial|r+'^owfc(A)|^trial) 
(Atrial I V'trial) 

which we minimize with respect to the parameters in "Atrial- Alternatively, we 
get a variational estimate by diagonalizing T+ Viowfc(A) in a truncated basis, 
where the trial wavefunction is a linear combination of the basis functions. 



3 Variational Results for the Deuteron and Triton 



Here we retrace the calculations of Ref. [2], starting with a study of the 
deuteron binding energy. As noted there, for weakly bound states we expect 
that a simple, generic ansatz should work increasingly well as the cutoff is 
lowered. Two such ansatze were considered. In the first one [11], the (un- 
normalized) and ^Di trial functions for the deuteron are (following the 
conventions of Ref. [12]) 

1 ak'^ 

= (F + f)(F + /.2) ' = WT7W+^' ' 

where 7, jj,, z/, and a are variational parameters. The underlying physics implies 
that /i and u should be roughly the inverse range of the interaction and 7 
should be close to (— MAri?^)^/^, where Mn is the mean neutron-proton mass 
and Ed —2.2246 MeV is the deuteron binding energy. The regulator in Viowfc 
implies that the corresponding exact deuteron wave function does not contain 
high-momentum components. Therefore, the two-body trial wave functions 
are multiplied by the same regulator f{k) for the relative momentum. 

We also adapted the form used for a high-accuracy representation of the 
deuteron wave function in Ref. [12], for which 

j/max Jmax 7^ 

Mk)^m)^^^^, Mk)^m)z^^, (11) 

where the rrij are fixed in a geometric progression: 

rrij = {-MNEdY/'^ + {j - l)mo , with mo = 0.9 fm"^ , (12) 

by treating the Cj and Dj coefficients as variational parameters for a given 
value of jmax- (The very accurate parameterization of the deuteron wave func- 
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Fig. 1. Deviation from of the best variational energy as a function of cutoff A 
for the wave function ansatze of Eqs. (10) and (11) with different numbers of terms. 
Smooth cutoff Viowfe results are given by the open symbols. 

tion for the Bonn potential in Ref. [12] has jma.x = 11, with some constraints 
on the Cj's and D^'s.) Since the variational coefficients appear linearly, we can 
simply diagonalize the Hamiltonian in the truncated basis of Eq. (11) to find 
the best estimate of the deuteron energy. 

The best variational energy for Eq. (10) as a function of a sharp cutoff is shown 
as the filled squares in Fig. 1. These estimates are not even bound for cutoffs 
above A pa 5 fm~^ (which includes the bare Argonne vis potential) but rapidly 
improve as the cutoff is lowered further, reaching a minimum deviation of 
less than 3keV around A 1.5 fm~^. Similar results are found for the ansatz 
of Eq. (11) with j'^ax = 3 (solid diamonds) and jmax — 4 (solid circles). As 
expected, lowering the sharp cutoff dramatically improves the effectiveness 
of the simple wave function ansatze, but the results unexpectedly worsen for 
cutoffs that are significantly larger than the naive limiting value set by the 
"binding momentum" of the deuteron, 0.25 fm~^. 

The variational estimates for the smooth cutoff of Eq. (3) are shown in Fig. 1 
as unfilled symbols for the corresponding ansatze. For the ansatz of Eq. (10), 
the smooth cutoff results are inferior at intermediate cutoffs but continue to 
improve monotonically at smaller cutoffs. For the other ansatz, the smooth- 
cutoff results are superior throughout and improve until reaching agreement 
with the exact result at the eV level. This is in accord with intuition for such 
a low-energy bound state, and it emphasizes that one works much too hard in 
calculating low-energy observables using conventional potentials that contain 
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Fig. 2. Momentum-space ^S*! deuteron wave function for two cutoffs and for the 
bare Argonne vis potential. 

strong high-momentum components. 

The reduced effectiveness of simple ansatze for the deuteron when using sharp 
cutoffs below 2 fm~^ can be understood by looking at the corresponding wave 
functions in momentum space. In Fig. 2, we show the exact deuteron wave 
functions in momentum space for both the smooth and sharp cutoffs at A = 
1.6 fm~^ and 2fm~^. The Viowfc wave functions remove the short-range/high- 
momentum behavior (e.g., the node just above 2fm~^ [13]) that is increasingly 
resolved at higher cutoffs, requiring finer and finer cancellations in the vari- 
ational integrals and a more correlated wave function. At A = 2fm~^, the 
sharp wave function is well behaved and an ansatz cutoff at the same momen- 
tum is adequate for a momentum-space variational calculation. (Even here, 
the abrupt cutoff creates problems in coordinate-space calculations, particu- 
lar for the ^Di component of the wave function.) But by A = 1.6 fm~^, one 
clearly sees a complicated endpoint behavior that will not be reproduced in 
simple variational trial functions. In contrast, the smooth regulator potential 
and corresponding wave functions do not have these artifacts. 

The extension from the deuteron to the triton in Ref. [2] was kept simple by 
using a truncated harmonic oscillator basis for a variational calculation with 
the two-body interaction only, which we repeat here for the smoothly regulated 
Viowfc potentials. The antisymmetric three- nucleon basis is generated from the 
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Fig. 3. Convergence of the binding energy of the triton from a direct diagonalization 
of the Hamiltonian in a harmonic-oscihator basis with fixed oscihator parameter 
^osc = 1.5 fm. Results for two cutoffs are shown. In each case the fihed symbols are 
for a sharp cutoff and the unfilled symbols are for a smooth cutoff. 

Jacobi coordinate oscillator states [14] 



{nlsjt-Uc]^j\)JT), 



(13) 



where (nlsjt) and (A/'£|j7'|) are the quantum numbers corresponding to the 
two relative Jacobi coordinates [e.g., k = i(pi — P2) and q = |(p3 — |(pi + 
P2))], and the basis is truncated according to the total number of oscillator 
quanta = {2n+l+2Af+C) < Nmax- Diagonalizing the intrinsic Hamiltonian 
in the truncated basis and minimizing with respect to the oscillator length 
parameter h provides a variational estimate to the true ground-state energy. 



The same pattern of the variational estimates for sharp cutoffs seen for the 
deuteron, namely improvement to a minimum followed by worsening estimates 
for very low A, was observed for the triton in Ref. [2]. However, an additional 
problem of convergence with the size of the harmonic oscillator basis, which is 
of greater practical importance, was seen as well over a wide range of cutoffs. 
For cutoffs around 2 fm~^, energies within about 200 keV of the accurate result 
from Fadeev calculations are reached with relatively small basis size, but then 
further convergence as the basis is increased is extremely slow (see the filled 
symbols in Fig. 3). In fact, extremely slow convergence beyond the 50 keV 
level is found for a wide range of cutoffs. One expects convergence at the keV 
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Fig. 4. Deviation from of the optimized variational energy as a function of cutoff 
A for a harmonic oscillator basis of fixed size N^^x = 40, comparing a sharp cutoff 
to an exponential regulator with n = 4. The oscillator parameter is optimized for 
each A. 

level for good variational calculations [3] . 

The poor convergence of the sharp cutoff Viowfc in an oscillator basis is anal- 
ogous to "Gibbs overshoot" phenomena found in Fourier series expansions of 
discontinuous functions. This connection is supported by noting the rapid con- 
vergence of the energy with the smooth cutoff Viowfc (e.g., the unfilled squares 
for A = 2.0 fm"^ in Fig. 3 are within 2 keV at A^max = 20 and within 1 keV 
at iVmax = 30 of the converged result). Similar dramatically improved conver- 
gence using the smooth cutoff Viowk is also found at other cutoffs (e.g., the 
unfilled circles in Fig. 3). 

Fadeev results with smooth cutoffs are not yet available to verify that they 
agree with the converged energies obtained here. However, we can look back 
at deuteron calculations using a harmonic oscillator basis to anticipate what 
we will find. Figure 4 shows the deviation of the variational estimate from the 
exact deuteron binding energy for a harmonic oscillator basis, comparing the 
two types of regulator. Each point was optimized with respect to the oscillator 
parameter in a basis of fixed size A^max = 40. As found for the other ansatze, 
the smooth cutoff improves steadily with decreasing A. We expect similar 
behavior for the triton. ^ 

^ We do not expect the addition of consistent three-body forces, which have been 
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4 Conclusions 



In summary, use of a smooth cutoff regulator preserves the improvement with 
lower cutoff and remedies the technical problems noted in Ref. [2] for varia- 
tional calculations of the dcutcron and triton with low-momentum potentials. 
While the deuteron and triton arc not stringent tests for heavier nuclei, these 
results coupled with the rapid convergence of the particle-particle channel ob- 
served in nuclear matter [1] imply that low-momentum potentials will be much 
more effective for few-body and many-body variational calculations than any 
conventional large-cutoff potential. Investigations of the effectiveness of run- 
ning a smooth cutoff lower for chiral EFT potentials, which are themselves 
low-momentum potentials compared to conventional potential models, will be 
explored elsewhere along with many other issues [4] . 

To take advantage of these observations, the variational calculations described 
recently in Ref. [3], which show the advantages of nonlocal interactions, should 
be particularly well suited. Based on our results here, we anticipate even more 
efficient variational results for low-momentum interactions with smooth cut- 
offs, with the added advantage of being able to vary the cutoff as a tool to 
optimize and probe the quality of the solution. Furthermore, wc can avoid the 
problem of constructing consistent, model-independent operators for conven- 
tional potentials by concurrently evolving to low momentum the potential and 
operators from chiral EFT. 

Since Hartree-Fock becomes a reasonable starting point for nuclear many-body 
calculations [1], the large arsenal of techniques developed for the Coulomb 
many-body problem becomes available and should be explored as well. In 
addition, the development of the smooth cutoff Fiowfc allows one to unam- 
biguously perform a density matrix expansion starting from low-momentum 
interactions to gain insight into the microscopic foundation of the nuclear en- 
ergy density functional, as the required coordinate space quantities are now 
well-defined with the smooth cutoff [15]. 
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shown to become more perturbative with lower cutoffs [9,1], to alter this assessment. 
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